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Abstract 

The acceleration of charged particles in the presence of a magnetic 
field and gravitational waves is under consideration. It is shown that 
the weak gravitational waves can cause the acceleration of low energy 
particles under appropriate conditions. Such conditions may be sat¬ 
isfied close to the source of the gravitational waves if the magnetized 
plasma is in a turbulent state. 
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1 Introduction 

The interaction of a charged particle with a gravitational wave (GW) has 
been studied by using a Hamiltonian formalism (P-i)- This approach is 
similar to the one used for the study of the interaction of ions with electro¬ 
static waves in magnetized plasmas ( 0 - 0 )- The existence of chaos in phase 
space has been shown as a consequence of the overlap of resonances. Re¬ 
cently, the acceleration of the charged particles due to the chaotic diffusion 
has been associated with bursts |H]. However, such a diffusion is obtained for 
relatively high energy particles and becomes efficiently present only for large 
amplitude values of the GW (a > 0.1). 

Furthermore, the existence of a parametric resonance in the general prob¬ 
lem of interaction between particles and waves is discussed in the classical 
book of Landau and Lifshitz pj. For the particular problem Kleidis et al. |1] 
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showed that a parametric resonance arises for a set of frequencies of the GW 
of the form i/ = 2/n, n = 1, 2,... and an estimation of its width is given. Also, 
the parametric resonance has been associated with the existence of chaotic 
motion. 

In the present paper we study the possible fast acceleration of low energy 
particles for small amplitudes of the GW by taking into account two special 
properties of the system. The first is the parametric resonance which, as we 
shall show, arises only for v = 2. The second property is the integrability of 
the system, and subsequently the absence of chaos, when the GW propagates 
in a parallel direction with that of the magnetic held. 

In section 2 we review the Hamiltonian formalism of the model which 
describes the interaction. In section 3, we study systematically the rise of a 
parametric resonance and in section 4 we prove the integrability of the system 
in the case of parallel propagation. It is shown in section 5 that when the 
above integrable case is in parametric resonance fast acceleration of particles 
occurs. The astrophysical implications are discussed in section 6. 

2 Review of the system. The Hamiltonian 
formalism 

The dynamics of a charged particle in a curved space-time and in the presence 
of a uniform magnetic held can be studied by considering the Hamiltonian 
(see dm) 

H{x'^,Pk) = eAi){pj-eAj) = z, j, A; = 0,..., 3, (1) 

where and pk denote the generalized coordinates and momenta of the 
particle, respectively. The metric dehnes the curved space-time and 
Ai = Ai[x^) are the components of the vector potential A of the magnetic 
held. It can be considered that the curvature of space-time is due to a grav¬ 
itational wave (GW). By considering a linearly polarized GW propagating 
in a direction k of angle 6 with respect to the axis 2 ;, a possible form of the 
tensor is the following HD cm 

a < 0, (2) 

where = {1, —1, —1, —1} is the Minkovski space-time and a is the ampli¬ 
tude of the gravitational wave. By using the metric tensor derived in and 
expanding its components up to 0{a) we obtain the nonzero components of 
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the tensor h'’^ in the following form 

cos^ e, ^22 = sin2 e, sin(20), 

^ (3) 

where 

'h = cos sin 9 + cos 9 — x^)) . 

The constant parameter z/ is associated with the freqnency oj of the GW. 
In the case of a nniform magnetic field B = Bq^z, which corresponds to 
A = {0,0, 0}, we may choose a system of nnits where the Larmor 

freqnency of the particles is G = eB^/mc = 1. In this case v is the ratio of 
the GW freqnency to the Larmor freqnency i.e. v = uj/VL. 

By snbstitnting in m the above specihc expressions and performing a 
canonical transformation to new canonical variables {(j),q,J,p) (for details 
see IS]) we get 

Lf = ^ (j2 — 2(/ + psin6*) J — {p^ + q'^)) + ~ cos^ 9) cos(z/0) = ^ 

(4) 

The variables p, q are associated with the gyro motion of the particles in the 
magnetic field, J = 1 — 7, where 7 = (1 — is the particle’s energy and 

V its velocity. The qnantity / = P 0 +P 3 /cos 9 is a constant of motion dne to 
the symmetries of the original system. 

3 The basic periodic solution, stability and 
parametric resonance 

For the Hamiltonian dH) we can easily obtain the periodic solntion 

p{t) = q{t) = J{t) = 0, 0(f) = —It + 00 (mod27r) (5) 

The solntion © corresponds to rest particles since it provides 7 = 1 (the 
energy of the rest mass). Becanse of the restriction iL = 1/2 we shonld set 
J = 1 in order onr system to possess the above special solntion. Thns, the 
Hamiltonian is now written in the form 

H = J + p J sin 0 + ^ (p2 -I- _ ia(g2 — (.Qg2 cos(z/0) = 0 (6) 

and the solntion © is periodic of period T = 27r. 

By considering the variations ■Ci,■^ 2 ,Cs)'Cr corresponding to the variables 
0, q, J and p, respectively, the variational eqnations for the periodic solntion 
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Q are 


= sin6'.^ 4 , .^2 = sin 0 ,^3 + (1 + a cos^ 6^ cos(i^t)).^ 4 , (7) 

is = 0, ^4 = (-1 + a cos(z/t ))^2 

Furthermore, let M be the monodromy matrix. Considering the unper¬ 
turbed system (i.e. a = 0) we obtain from the monodromy matrix M the 
stability index 

k = 2 — trace(M) = 2 — 4 cos j (8) 

The condition —2 < k < 2, ensures the linear stability of the periodic orbits 
for all values u except for u = 2jn,n = 1,2,.... These values are candidates 
for the generation of a parametric resonance. 

In the perturbed system we examine the stability by following the method¬ 
ology described by Yakubovich and Starzhinskii PI. p.357. The variational 
equations (jH)) are of the form 

e'=(Ko + aDi(f))e (9) 

and the fundamental 4x4 matrix of solutions is 


S(f;a) = , K(a) = Kq + aKi + + ... (10) 


where F(t; a) is a T-periodic matrix, Kq is the matrix of coefficients of the 
linear system (jH]) and Kj, i = 1,2,... are constant matrices calculated from 
the variational equations (jH)) on the periodic orbit. After some calculations 
we hnd K* = 0,Vi ^ 0 , 2 . The nonzero components of the matrix K 2 = 
{^ 2 (h j)} are 


772 ( 2 , 3 ) 


1 c(362 + 1/2-1) 

2 z/4-51/2 + 4 ’ 


^ 2 ( 2 , 4 ) 


1 6(1 + 6 ) 

2 1/2 - 4 ’ 


^ 2 ( 4 , 2 ) 


1 (1 + 6 ) 

2 j/2 - 4 ’ 


where b = cos^ 6 and c = 6 sin 6 '. The stability of the periodic solution is 
given by the eigenvalues of the matrix K(a) dehned in ©, which are 


Ai — A2 — 0, A3 


—A 4 = i 


(. a\l + bf \ 
V 4(^/2 - 4 ) ) 


+ 0{a^) 


( 11 ) 


Conclusively for every value of a,d and v ^ 2 the eigenvalues A 3,4 are purely 
imaginary and the periodic orbit that corresponds to the rest particles is 
stable. Therefore, for arbitrarily small amplitude values a, instability (para¬ 
metric resonance) may occur only at i/ = 2 . At this value the above method 
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Figure 1: The region of the parametric resonance in the parameter space 
1 / — a for some typical values of the propagation angle 9. 


fails to classify the periodic orbit and we proceed to a numerical investiga¬ 
tion. Namely, we obtain the monodromy matrix by integrating numerically 
the variational equations beside the canonical equations of the Hamiltonian 
(|F)jl for various values of the amplitude a of the gravitational wave and of the 
propagation angle 9. 

As it is shown in hgure 1, the periodic orbit becomes unstable in a V- 
shaped region of parameters centered at z/ = 2. We observe that the borders 




Figure 2: Poincare sections of the phase space around the periodic orbit 
(g = p = 0) for a = 0.01, 9 = 30° and (a) u = A, (h) u = 2, (c) u = 1. 
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of the parametric resonance can be approximated by the fitting formula 

a = (1 + sin^ 6)Ai>, Au = \i> — 2\ 

and the instability region broadens as 6* —0, i.e. towards the parallel propa¬ 
gation. In hgure 2 we present the Poincare sections of © on the plane q—p 
for 0 = 0 (mod(27r/i/)). The periodic orbit correspond to the hxed point 

of the sections. We see that at the parametric resonance a.t u = 2 a pitch- 
fork bifurcation of the periodic orbit © occurs. Two stable periodic orbits 
bifurcate, which correspond to 7 7^ 0, while the periodic orbit © becomes 
unstable and the associated stable and unstable manifolds form a separatrix 
curve of figure-eight. 

The Poincare section of system © can be defined also in the plane of 
variables 0 — 7 for q = 0, H = 0. Then the topology of the section in Fig.2b 
is represented in the section 0 — 7 as it is shown in Fig.3a. The unstable 
periodic orbit B. indicated by Pq, appears at 0 = ±7r/2 and its separatrix 
curve is located in the enegry domain 0 < 7 < 75. For any initial conditions 
the particle’s energy 7 shows regular bounded oscillation of librational type 
around Pi or rotational type above the separatrix S. The stable periodic 
orbit, which is generated by the pitchfork bifurcation and indicated as Pi, 
corresponds to energy value 71. As it is illustrated in Fig.3b, the value 71 
depends on both parameters a and 9. It holds always 71 < 75 and also we 
can obtain by fitting that 71 ~ 6~‘^ for 0 —> 0. Consequently, as 0 —0, 
the energy oscillations increase rapidly their amplitude and at 6* = 0° all 
particles are ejected to infinite energy independently of the amplitude of the 
GW. This special case is studied in the following section. 

4 The case of parallel propagation 

We consider 6 = 0°, i.e. the GW propagates in parallel direction with the 


magnetic field. Then the Hamiltonian © takes the form 

H = J + ^{p‘^ + q'^) + - q^) cos(z/0) = 0 (12) 

The corresponding canonical equations are: 

0 = 1 (13) 

j = (^4) 

q = p(l-|-a cos 120) (15) 

p = q'(—1-|-a cos 120) (16) 
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Figure 3: (a) Poincare section presented in the 0 — 7 plane for u = 2, a = 0.01 
and 9 = 30°. (b) The energy 71 of the bifurcated stable periodic solution as 
a function of the propagation angle 6. As 6 ^ —0 we obtain that 71 ~ 6~‘^. 

From equation m we obtain directly that 0 = t + 0 o and subsequently 
equations m and (uni) consist a linear non-autonomous system. Providing 
a solution p = pit), q = g(t) the system turns to be integrable for any values 
of the parameters a and z/. The Poincare section in Fig.4a, indicates the 
regular oscillations of the particle’s energy for all initial conditions. Such a 
situation changes at the parametric resonance for v = 2 and the dynamics is 
depicted by the Poincare section in Fig.4b. We obtain that the stable and 
the unstable manifolds which originate from the unstable periodic orbit © 
are located at 0 = (2fc —l) 7 r /4 and (j) = (2A; + l) 7 r /4 [k ^ Z), respectively, and 
extend from 7 = 1 up to inhnity. The dynamics becomes purely hyperbolic. 

The linear system of eqs. m and (uni) is not directly solvable. By 
differentiating and introducing the variable transformation u = p + q and 
V = p — q, we derive the separable system of second order equations 

ii + (1 + az/sin(z/t))u = cos^(z/t)M (17) 

h + (1 — az/sin(z/t))u = cos^(z/t)u. (18) 

If we introduce the new time variable r = 7 r /4 — z/t/2 and omit the terms 

Oia?) the above equations take the form 

X" + (^4 ± — cos(2r)^ X = 0, (19) 

where x stands for either u 01 v, the minus sign refers to v and the plus sign 
to u and prime denotes differentiation with respect to r. Equation (P|) is 
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Figure 4: Poincare sections in the case of parallel propagation {6 = 0°) for 
(a) z/ = 5, a = 0.1 and (b) z/ = 2, a = 0.001. The up and down directed 
arrows indicate the unstable and stable manifolds, respectively. 

a Mathieu type equation with solutions of the form x(t) = e^'^tc(r), where 
w { t ) is a periodic function jHj. In the parameter space z/ — a stable and 
unstable regions are defined where Re[X\ = 0 and Re[X] ^ 0 respectively. 
For the actual system m the unstable region is the V-shaped region shown 
in Fig. 1 for 6 = 0°. In Fig.5a we present an example of the oscillations of 
the variable q which correspond to the trajectory along the unstable manifold 
for u = 2. Such a variation is associated with an increase of the particle’s 
energy 7 = 1 — J. From m it follows that 

7 = 1 + + 9") + 0(a) = 1 + i(«" + v^j + 0 ( a ) (20) 

The increase of 7 is exponential (see Fig. 4b) and the numerical study suggest 
that for f 3> 1 , the associated exponent coincides with the amplitude of the 
GW, i.e. 

7 ~ for t 3> 1 ( 21 ) 

5 Passage through the parametric resonance 

In previous studies, mentioned in the introduction, the acceleration of the 
particles is a result of the chaotic diffusion. However, wide chaotic regions 
are formed only for high energy levels and in the presence of very strong 
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Figure 5: (a) The increasing amplitude q = q{t) of the gyro motion of a 
particle in the parametric resonance (b) the exponential increase of energy 
7 = 7 (f) for the evolution along the unstable manifold. The corresponding 
parameters are u = 2, a = 0.01 and 7 ( 0 ) = 2. 

gravitational waves (a > 0.01) [S]. Especially, for small values of the propa¬ 
gation angle 6^0 chaotic regions appear in phase space only as very thin 
bounded zones around islands formed according to the Poncare-Birkhoff the¬ 
orem. Consequently, acceleration of low energy particles under the perturba¬ 
tion caused by weak GWs can occur only through the parametric resonance 
(z/ = 2) at the integrable case {9 = 0°). 

In a more realistic situation, the direction and the magnitude of the mag¬ 
netic held met by the particles is not constant. According to the results 
obtained in the previous sections, the charged particles will accelerate when 
the magnetic held becomes almost parallel to the direction of the gravita¬ 
tional wave and have a magnitude such that v k, 2. A simple simulation that 
illustrates this situation is described below. 

We consider that the propagation angle between the GW, with constant 
u = 2, and the direction of the magnetic held varies in time linearly as 

e{t) = eo + Xt ( 22 ) 

The time variable t is the coordinate time in the original form (HD of the 
system. The slope parameter A controls the time span spent by the particle 
in an interval (^ 1 ,^ 2 ) and, particularly, the critical time interval At* where 
0 0. A typical example of the evolution of a particle under the above 

condition is shown in Fig. 6 a. 

We observe that the particle, which starts with low energy 7 = 70 (par¬ 
ticularly 7 o = 3) and 60 = —5°, obeys small oscillations for approximately 
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Figure 6 : (a) The evolution of 7 along a trajectory in the system with vary¬ 
ing propagation angle, v = 2 and a = 0.001. The energy jump is clearly 
seen when 0 0°. (b) The magnitude of the energy jump Ay against the 

amplitude of the GW. 

9 < —1°. When the system approaches the integrable hyperbolic dynamics 
at 6 * = 0 °, the particle accelerates and a signihcant jump in particle’s energy 
occurs. Leaving the hyperbolic domain {9 > 1) the particle is captured in a 
regular phase space domain with corresponding energy 7 = 71 signihcantly 
greater than 70 (particularly, 71 Ri 20). The total increase in energy. Ay, 
depends on the critical time interval At* and the 6 '-width, A^*, where the 
hyperbolic structure is effective. For the particular simple model it holds that 
At* = A9*/X. Certainly, the most signihcant parameter for the acceleration 
is the amplitude a of the GW which is associated with the exponential rate 
of the acceleration along the unstable manifolds in phase space according to 
the estimation m- In Fig. 6 b we present the relation between a and Ay 
derived by the numerical simulation for some particular values of the param¬ 
eter A. We observe that, if the critical time At* is sufficiently large (or A is 
small), the energy jump is notable even for very small values of a. The same 
qualitative results are obtained when we keep constant the propagation angle 
9 = 0° and change the parameter u passing through the value v = 2. 


10 




















6 Astrophysical implications and conclusions 

We have shown that gravitational waves with frequency lo and amplitude a 
can accelerate charged particles to very high energies when two conditions 
are met (a) the gravitational wave propagates along the external magnetic 
held [d ~ 0°) and (b) the ratio of the gravitational wave frequency to gyro 
frequency is about 2 (i/ 2). The model used has a number of important 

simplihcations: (1) We assume that the magnetic held is constant (2) The 
energy absorption by the particles is small compared to the energy curried 
by the gravitational wave and (3) we ignore the collective plasma ehects 
associated with the particle acceleration. 

Retaining the above simplihcations we may proceed to estimate the en¬ 
ergy gained by the particles in a realistic astrophysical system. We assume 
that a volume of length L lO^^cm is close to a source of gravitational waves. 
The volume is hlled with plasma with density n ~ 10^^particles/cm“^ and 
the magnetic held is in a fully turbulent state. The gravitational wave is 
crossing the volume in few seconds and every time it travels along the am¬ 
bient held {9 K. 0°) with almost zero strength its energy increases. The zero 
strength magnetic held, which is associated to null surfaces, is an important 
requirement since the gravitational wave has low frequency (of few KHz) 
and acceleration is possible only for v k. 2). The energy increase is fast and 
can be substantial (Ay 10) if the amplitude of the gravitational wave is 
a ~ 10“'^ — 10“®. Therefore the particles dihuse in energy space not by small 
random steps, as it was the case in the stochastic interaction studied earlier 
(0- i), but through large random energy jumps. 

It is beyond the scope of this article to present a detailed analysis of 
the above process but we can present a rough estimate of the energy trans¬ 
ferred from the gravitational wave to the plasma. Let us assume that a 
particle travel a distance Lp > L (its trajectory is not a straight line) be¬ 
fore exiting the interaction volume. Let us also assume that only in small 
parts of its trajectory {Lint ~ a few kilometers) the particle is able to meet 
the conditions needed for its acceleration. For a small portion of protons, 
e.g Uacc = '^totai/100), wliich participate in the interaction, these regions 
can be a fraction of N^ax = Lp/Lint ~ 10^ and the energy gain will be 
Ekin ~ 10 ®(A 7 — l)mpC^ lO^ergfs (or lO^^eK), where rup is the mass of 
proton. The number of particles participating in this interaction can be as 
many as Ntotai ~ naccL^ ~ 10"^®. The total energy lost by the gravitational 
wave is, approximately, 10^® — 10“^® ergs, which constitute a small fraction of 
its total energy 10^® ergs (see na). 

According to the model, the acceleration conditions {6 ^ 0°, u ^ 2) can be 
met sporadically along the trajectory of a charged particle in the vicinity of a 
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gravitational wave source when the magnetized plasma is in a turbulent state. 
A few particles can then reach very high energies in a few seconds absorbing 
only a very small fraction of the gravitational wave. It remains to be seen, 
by performing further simulations, if the process described above for particle 
acceleration can be significant in a realistic astrophysical environment. In 
this case the electromagnetic emission from the accelerated particles can be 
a precursor of the gravitational wave for easier detection. 
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